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Abstract. 

We summarise results for the predicted density and peculiar velocity fields 
from the PSCz survey, consisting of redshifts for 15,000 IRAS galaxies covering 
84% of the sky to a depth of 25, 000 km s^^. We have used a generalisation 
of the Path Interchange Zeldovich Approximation technique to determine the 
velocity field; the most remarkable feature being a coherent large-scale flow 
along the baseline connecting the Local Supercluster, Centaurus and the Shapley 
Concentration. Comparison of the predicted and observed bulk flows gives a 
value of /3 = 0.50 ±0.1. 

We re-examine the PSCz dipole, with improved redshift completeness at 
low-latitudes and using PIZA to estimate real-space distances. We find the 
dipole to be stable between 80 and 180 ft.~"'^Mpc, although there appears to be a 
significant contribution to the dipole around 200/i^^Mpc. The overall misalign- 
ment with the CMB dipole remains at 20°. The implied value of /3 depends on 
the exact treatment; we derive values in the range (3 = 0.40 — 0.55 with statistical 
error ±0.1. 

We also present the density field and a preliminary dipole analysis from the 
almost completed Behind The Plane survey, which extends the PSCz to cover 
93% of the sky including the proposed core of the Great Attractor. We find 
a density peak at (325, — 5, 3500kms~^), about half as massive as Centaurus 
or Pavo-Indus, and forming part of a continuous filament linking them across 
the Plane. We also find evidence for a much larger 'Greater Attractor' directly 
behind the GA, at a distance of 125 h~^Mpc, and more massive than the Shapley 
concentration. 

At large distances the dipole direction is in much improved agreement with 
the CMB: at 250 - 300/i~iMpc, the misahgnment is only 5 - 10°, and this is 
consistent with the shot noise errors. The Ja-weighted dipole, which in principal 
is expected to agree better with the CMB, gives a misalignment of 13°. The 
dipole direction is somewhat dependent on the corrections made for the redshift 
incompleteness, and may also be affected by unquantified incompleteness in the 
BTP survey close to the Galactic Centre. The dipole amplitude implies a value 
for the parameter (3 of 0.44 ± 0.1. Values much larger than /3 — 0.5, while not 
formally ruled out, cause very unlikely-looking rocket effects. 
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Figure 1. PSCz and BTP galaxy catalogues and masks in galactic coordinates. 



1. Introduction 

Redshift surveys of galaxies can be used to infer the dynamics of the Local Universe. 
Comparison with observed peculiar velocities allows determination of the mass associ- 
ated with galaxies (usually parameterised by the the quantity /3 ~ fl^.d/b), and also 
reconstruction of the initial conditions. Because gravity is extensive, accurate predic- 
tions formally require uniform surveys with infinite depth and complete sky coverage; 
failing that, we need a depth much larger than the region of interest, and with coverage 
gaps small enough to interpolate over. IRAS surveyed 96% of the sky with negligible 
extinction, tracing the large-scale structure to a depth of at least 20, 000 km s~^. For 
these reasons, IRAS-based surveys have dominated this field for the last decade. The 
IRAS dipole has remained astonishingly stable in direction over this period, despite 
the ever-increasing depth and sky coverage of the follow-up redshift surveys. However, 
the amplitude has increased continuously, with extra contributions uncovered at larger 
distances and lower latitudes. 

The Point Source Calatog Redshift Survey (PSCz) consists of essentially all galax- 
ies in the PSC, down to its nominal flux limit of 0.6 Jy, across the 84% of the sky where 
extinction was estimated to be low enough {As < 2™) to allow reasonably complete 
optical identificati ons and spectroscopy with 2-mctrc class telescopes. T he PSCz data is 



now available from http: / /www-astro. physics.ox.ac.uk/~wis/pscz.html ; also there avail- 



able is the full, labelled 3D density distribution, using the Fourier interpolation method 
presented in Saunders and Ballinger (2000). The catalogue and main science papers to 
date are summarised in Saunders et al. (2000a). 

The PSC itself is useable much deeper into the Plane. The Behind The Plane 
Survey extends the PSCz to cover the entire 93% of the sky with useable PSC data, 
and is now almost completed. The remaining size of the Zone of Avoidance is just 4% 
of the sky, and this is hopefully small enough to interpolate over with confidence. The 
sky distribution for both surveys is shown in Figure 1. 
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2. The PSCz velocity field and the Path Interchange Zeldovich Approxi- 
mation 

A linear theory reconstruction of the velocity field and comparison with the Mark III 
dataset was been performed by Branchini et al. (1999). However, the heavy smoothing 
needed to make linear theory applicable seriously compromises the detail with which 
such comparisons can be done. 

The PIZA algorithm, originally presented by Croft and Gaztaiiaga (1997), is a 
quasi- linear method for reconstructing the velocity field from galaxy surveys, allowing 
the full resolution of the data to be kept. The basic idea is wonderful and very sim- 
ple: to Zeldovich order, velocities and trajectories depend in a fixed way on the scale 
factor; so minimising the action is equivalent to minimising the rms trajectories of the 
galaxies. Given a set of uniform initial positions and the current galaxy positions, the 
least-action solution is found by setting up arbitrary trajectories and swapping pairs to 
reduce the action. The formalism can be readily adapted to minimise (for a given /3) 
the real-space trajectories from redshift-space information. With more complexity, the 
changes in selection function for individual galaxies caused by peculiar velocities can 
also be included, though we have found doing this leads to instabilities without any 
improvement in accuracy, and this is not done here. This and edge effects lead to a 
systematic velocity underestimation of ^ 20% in comparison with simulations. Further 
details are presented in Valentine, Saunders and Taylor (2000). The resulting veloc- 
ity field in the Supergalactic Plane is shown in Figure 2. The flow field is dominated 
by infall into the Supergalactic Plane, and coherent flow along it involving the Local 
Supercluster, Centaurus and the Shapley Concentration. 

3. Comparison with observed peculicir velocities 

The velocity field predicted from PIZA can be compared with observations. The derived 

trajectory of the Local Group is robust, with weak dependence on how the mask is filled 
or the input /3. Comparison of this trajectory with the CMB dipole gives (3 = 0.53 ±0.1. 
A comparison of the bulk flow as a function of radius with the Marklll dataset implies 
a value for /3 of 0.5 ±0.15. 

We have also undertaken a new comparison with the SFI dataset of peculiar veloci- 
ties of Spiral Galaxies (Giovanelli et al. 1998). We have calculated the average peculiar 
velocity with respect to the Local Group, of galaxies in redshift slices, so as to be as 
closely as possible comparable to Giovanelli et al. . We find reasonable agreement in am- 
plitude, though there are significant directional differences, for (3 = 0.55 ±0.1 (Table 1). 
We agree with Giovanelli et al. that most of the dipole is generated within 5000 km s~^. 
However, we do find substantial contribution to the dipole beyond this distance (see be- 
low), even though our derived velocity dipole at this distance is in excellent agreement 
with the CMB. The sum derived by Giovanelli et al. is a galaxy-weighted vector sum, 
rather than volume-weighted. Therefore, a good agreement between velocity vector and 
CMB dipole does not necessarily imply that a cosmic rest frame has been reached. 



SheU 


VsFI 


{Lb)sFi 


VpiZA 


{l,b)pizA 


0-2000 


270 ± 80 


(245, 49) ± 19 


338 


(237,41) 


1500-3500 


410 ±69 


(255, 21) ± 12 


404 


(265,51) 


2500-4500 


620 ± 76 


(255, 15) ±11 


471 


(271,45) 


3500-5500 


585 ± 92 


(265, 19) ±13 


550 


(285,36) 


4500-6500 


544 ± 98 


(270, 16) ±15 


508 


(284,31) 



Table 1. SFI vs PSCz peculiar velocities with respect to shells of galaxies. PSCz 
reconstruction assumes (3 = 0.5. 
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Figure 2. Reconstructed velocity field in the Supergalactic Plane. The Great Attractor 
region is at 10 o'clock. 

4. PSCz dipole 

In the course of the BTP survey, we also took several hundred new reshifts for low- 
latitude PSCz galaxies, giving an overall redshift incompleteness for bright galaxies of 
70 galaxies (0.5%). A further 230 galaxies (1.5%) without redshifts are assumed from 
their optical and/or infrared appearance to be beyond 30, 000 km s^^. In the original, 
publicly available survey, there was known incompleteness at lower redshifts towards 
the Galactic anti-centre as discussed by Saunders et al. (2000a) and whose effect on 
the dipole is discussed in Rowan-Robinson et al. (2000). It is appropriate then to 
recalculate the dipole analysed by Rowan-Robinson et al. . We have corrected for 
peculiar velocities using PIZA as above, with (3 — 0.5 for self-consistency between 
gravity and velocity dipoles. The mask is filled in via the Fourier interpolation method 
of Saunders and Ballinger (2000). The revised dipole is presented in Figure 3. The 
dipole is remarkably stable between 125 — 180 /i~^Mpc, However, it shows a significant 
shift at 180 — 220 h^^Mpc in the X and Y components. This shift is seen in the analysis 
of Rowan-Robinson et al. , and is too large to be attributed to the small remaining 
redshift incompleteness. It remains when the interpolated galaxies in the Plane are 
replaced by randomly distributed galaxies at the mean density, so it is not due to the 
interpolation either. It does not appear to be associated with any single particular 
structure. 

In linear theory, the contribution to the velocity dipole per galaxy is 



V(r) = 



4tt r^^{r) 



(1) 



where ip is the selection function. In Figure 4, we plot the amplitude of this velocity as 
a function of distance for the PSCz. The dipole contribution per galaxy is remarkably 
constant in the range 20 — 100 ft.~^Mpc. Nearby, the contribution is large because of the 
1/r^ term, but at least the numbers of galaxies involved are small and the structures 
they reside in are real and well-sampled. At large distances, the contribution blows up 
because tp falls off much faster than 1/r^, and the dipole becomes increasingly dominated 
by shot noise (e.g. Taylor and Valentine 1999). Strauss et al. (1992) derived an optimal 
weighting scheme, in the sense of producing the minimum variance ensemble difference 
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PSCZ cumulative dipole amplitude 

PSCz cumulative dipole direction 
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Figure 3(a). X,Y,Z and total velocity components of the real-space PSCz dipole, with 
redshifts corrected for peculiar motions using the PIZA reconstructions, with an as- 
sumed value of [3 — 0.5. The dashed lines show the ,J;^-wcightcd components. 
Figure 3(b). Angular components in Galactic coordinates; showing direction at 100, 200 
and 300/i~^Mpc (filled dots) with their one-sigma shot noise (dotted circles), also the 
Js-weighted direction (cross) and associated one-sigma shot noise (dashed circle), and 
the CMB dipole (star) and circles showing 5° and 10° deviation. 



between true and estimated gravity vectors. At large distances, this amounts to giving 
each galaxy a ' J3 weighting' i.e. 

47r ^ 1 + JsirWr) r^^(r) ^ ' 

gals 

Where Jz{r) = Air r'^^{r')dr' is the integral over the correlation function. Jsir) 
is not well-determined on large scales, but the precise value used is not critical as 
any error only means the weighting is quite optimal, while leaving the dipole estimate 
unbiased in an ensemble sense. Based on the PSCz power spectrum (Sutherland et al. 
1999), we have simply assigned Jsir) a constant value J3 = 20, 000 ft-^^Mpc for all r. 
The resulting velocity dipole per galaxy is shown in Figure 4; there is now a sharp 
attenuation above 200/i~^Mpc. 

Based on Lincweavcr et al. (1996) and Yahil, Sandage and Tamman (1977), we have 
assumed a value for the LG motion with respect to the CMB of (276, 30, 627kms^^). 
The misalignment of CMB and IRAS dipoles remains 15 — 25°, depending on chosen 
upper cutoff radius for the gravity sum, precise treatment of the mask etc. The inferred 
value of /3 is /3 = 0.54 ± 0.10 from the overall amplitude; however, the Z-component, 
which is very stable and also less sensitive to the inferred structure behind the Plane, 
gives P = 0.40 ± 0.07. 

The dipole and peculiar velocity results are largely independent, since the latter is 
sensitive only to tidal effects from large scales. The completely independent analysis of 
Tadros et al. (1999), using the statistical distortion of the clustering in redshift space, 
also gives a consistent, low-/3 result /3 = 0.47 ± 0.16. 



5. The Behind The Plane survey 

The 84% sky coverage of the PSCz survey is effectively limited by the need to get, 
for every galaxy, an optical identification from sky survey plates. The fractional incom- 
pleteness in sky coverage translates directly into a much larger uncertainty in predicting 
the gravity dipole on the Local Group: suppose the fractional masked area is /, and also 
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Infall per galaxy vs distance for PSCz 
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Figure 4. Linear theory velocity generated at the Local Group for a galaxy in the PSC 
for a galaxy of given distance. The dashed line shows the effect of the Js-weighting. 



that this is large and compact cmough that wc have essentially no information about 
the distribution behind the masli. Tlicn we can think of the sky as consisting of 1/ f 
independent contributions to the dipole, one of which is unknown, giving a fractional 
uncertainty in the measured dipole vector of \J f /{I — /). Thus even 10% missing sky 
coverage can lead to a error of 1/3 in the dipole. The bulk of the dipole is generated 
around 40 ft,~^Mpc, at which distance the correlation length ro subtends just 5°, so gaps 
much larger than this arc too large to interpolate usefully and lead to large errors in 
the dipole. 

The IRAS PSC data itself is reliable to much lower latitudes, although genuine 
galaxies are outnumbered by Galactic sources with similar IRAS properties. Previous 
attempts to go further into the Plane have either been restricted to the Arccibo decli- 
nation range, or have relied on optical identifications from Sky Survey Plates. Because 
the extinction may be several magnitudes or more, they have inevitably suffered from 
progressive and unquantifiable incompleteness as a function of latitude. Hence in 1994 
we embarked on a program, parallel with the PSCz survey, to systematically identify 
and get redshifts for low latitude IRAS galaxies wherever the PSC data allowed - es- 
sentially limited only by IRAS coverage and confusion, to 93% of the sky. A detailed 
description is given in the associated poster paper (Saunders et al. 2000b). The iden- 
tification program is complete; and we have found over 200 galaxies which were not 
included in the 1.2Jy survey but should have been. The BTP sky distribution is shown 
along with the PSCz in Figure 1. At low latitudes towards the Galactic centre, the 
PSC is not complete to our nominal flux limit of 0.6 Jy, and the density of identified 
BTP galaxies reflect this. This incompleteness is interpreted, quantified and corrected 
for in the poster paper. Notwithstanding this, the huge overdensity of galaxies in the 
Great Attractor region (325, ±5) is obvious; as is the Puppis supercluster (250, 0) and 
the extension from Perseus (150,-10) and Pegasus (95,-10) into the Plane. 

The overall surface density of IRAS galaxies in the BTP area, when corrected for 
known incompleteness, is 78% of the PSCz average. Unquantified incompletenesses 
could conceivably push this up to the average value, but clearly we can rule out large 
overall overdensity behind the Plane. We have used the Fourier interpolation technique 
to find the projected surface density of the PSCz+BTP galaxy distribution; this is 
shown in Figure 5. In addition to the above features, we see the Ophiucus Spercluster 
clearly at (0,10) and the local void behind the Galactic Centre. 
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Figure 5. Surface density of PSCz+BTP galaxies in galactic coordinates. 
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Figure 6. Weighted N(z) distribution for the BTP survey. Open circles show the 

effect of assigning a nearest-neighbour redshift to galaxies without measured redshifts 
themselves. The line is the prediction from the high-latitude PSCz. 




\ai= 3571 km/s, p = O.Ol^^^^H 100 

Figure 7. Galaxy density field at a redshift of 3571 km s"^ in galactic coordinates, 
showing the peak of the Great Attractor at (325,-5). 

6. 3D distribution 

The spectroscopy is now completed in the south, and useably completed in the north. 
In addition to BTP galaxies, we also systematicall pursued faint, low-latitude PSCz 
galaxies where there was any likelihood that these were within 30, 000 km s"\ There 
are c.370 very faint galaxies in the combined PSCz-|-BTP with unknown redshift, but 
these are likely beyond 30, 000 kms~^. There are a further 240 for which we would like a 
redshift but do not currently have one, though for half of these we have unreduced data. 
As a temporary fix, for these latter 240 galaxies we have simply assigned the redshift 
of a nearby galaxy with redshift. The N(z) distribution, both with and without this 
fix, is shown in Figure 6, together with the prediction from the PSCz survey. There is 
some evidence that we are missing galaxies at most or all redshifts, but also that there 
is no cutoff in the redshift completeness with distance till at least 25, 000 km s~^ - much 
better than our original target of 15, OOOkms"^. 

We can then make a 3-dimc nsional map of the galaxy distrib ution. The full, an- 



imated 3D map is available from http://www.roe.ac.uk/willwww. Figure 7 shows the 



peak we find in the GA region at (325, —5, 3500 km s ^). The overdensity is quite mod- 
est compared with other supcrclusters in the region, and the peak overdensity is closer 
than the velocity of the Norma cluster at 4844 km s~^ (Woudt et al. 2000). However, 
we find the GA to be very extended in redshift space, from 3000 - 5500kms"\ More 
surprising is Figure 8, showing an apparent 'Greater Attractor', directly behind the 
GA at (326, —3, 12500 km s^^). Because of the incompleteness weighting, and because 
the overdensity is right on the edge on the mask, the finding represents only a dozen 
measured redshifts and is at this stage preliminary. 



7. The PSCz+BTP Dipole 

The combined PSCz-|-BTP sample is now complete enough to attempt a revised dipole. 
One of the main motivations for this is that current low 51 models predict small mis- 
alignments between gravity and velocity vectors, especially as the Hubble flow in the 
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Reconstructed PSCz+BTP density field 




Figure 8. Galaxy density field at a redshift of 12313 km s ^ in galactic coordinates, 
showing the 'Greater Attractor' at (326,-3). 



immediate vicinity of the Local Group (where non-linear effects are most likely to arise) 
is so cold: a misalignment as large as the 20 — 25° between IRAS and CMB dipoles would 
present problems for current models. As above, we have assigned a nearest neighbour 
redshift to 240 galaxies. We have used our interpolated 3D maps to generate random 
'galaxies' to add to the catalogue when determining the dipolc. However, because the 
missing region is so much smaller, the results are now much less sensitive to this inter- 
polation. 

To get real-space positions, we used the PIZA reconstruction for the PSCz galaxies 
as above, with an assumed P = 0.5. For the BTP sources we just used the redshift-space 
position. The resulting dipole is shown in Figure 9. Overall, the dipolc looks much the 
same as for the PSCz alone. However, the misalignment with the CMB is now very much 
smaller, being 10° at 250 h~^Mpc and 5° at 300 h~^Mpc, with the shot noise error being 
10°. The misalignment is 13° for the Js-weighted sum. The directional agreement also 
appears to be robust: if the galaxies with unknown redshift are dropped altogether, the 
dipolc only changes by 5° in direction and 7% in amplitude, while if random sources 
are put uniformly in the mask at the average density instead of interpolated, the dipole 
only alters by a few degrees in direction and 10% in amplitude. 

The shift at 180 — 220 h~^Mpc seen in the PSCz remains. At this stage, it remains 
possible that variable redshift completeness is responsible. Also, at this distance, very 
slight systematic effects in the PSC itself could have large effects on the dipole, but 
these would not be expected to be restricted to a particular velocity range as seen in 
Figure 8. Beyond that, there is a drift in the X component, but this too is sensitive 
to variable redshift incompleteness around the Galactic Plane and is not secure at this 
stage. 

Peacock (1992) noted that apparent convergeance in the dipole may be illusory. 
The apparent convergeance of the IRAS dipole between 80 — 180 h~^Mpc, in a direction 
only 20° away from the CMB, was seen in the QDOT, 1.2Jy and PSCz surveys, and 
was in each case interpreted as the 'true' value for the determination of /?, with no 
significant contribution from larger distances. We now sec that there may indeed be a 
significant contribution from beyond 180/i^^Mpc. Of course the same argument can in 
principal be applied this work, in that we have not ruled out contributions from beyond 
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Figure 9(a). X,Y,Z and total velocity components of the real-space PSCz+BTP dipole, 
with PSCz redshifts corrected for peculiar motions using the PIZA reconstructions, 
with an assumed value of /9 = 0.5. The dashed lines show the Js-weighted components. 
Figure 9(b). Angular components in Galactic coordinates; showing direction at 100, 200 
and 300/i^^Mpc (filled dots) with their one-sigma shot noise (dotted circles), also the 
Js-weighted direction (cross) and associated one-sigma shot noise (dashed circle), and 
the CMB dipole (star) and circles showing 5° and 10° deviation. 

300/i~iMpc. However, we can hope that the excellent directional agreement we find, 
combined with the implausible size of structures required, make this unlikely. 
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